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ABSTRACT 

We explore the origin of FUSE and HST STIS far UV spectra of the dwarf 
nova, EY Cyg, during its quiescence using combined high gravity photosphere 
and accretion disk models as well as model accretion belts. The best-fitting 
single temperature white dwarf model to the FUSE plus HST STIS spectrum of 
EY Cygni has Teff = 24, OOOK, log g = 9.0, with an Si abundance of 0.1 x solar 
and C abundance of 0.2 x solar but the distance is only 301 pc. The best-fitting 
composite model consists of white dwarf with Tg// = 22, OOOK, log g — 9, plus 
an accretion belt with T^ej* = 36, OOOK covering 27% of the white dwarf surface 
with Vheitsini — 2000 km/s. The accretion belt contributes 63% of the FUV light 
and the cooler white dwarf latitudes contribute 37%. This fit yields a distance 
of 351 pc which is within 100 pc of our adopted distance of 450 pc. EY Cyg has 
very weak C IV emission and very strong N v emission, which is atypical of the 
majority of dwarf novae in quiescence. We also conducted a morphological study 
of the surroundings of EY Cyg using direct imaging in narrow nebular filters from 
ground-based telescopes. We report the possible detection of nebular material 
associated with EY Cygni. Possible origins of the apparently large Nv/Civ 
emission ratio are discussed in the context of nova explosions, contamination of 
the secondary star and accretion of nova abundance-enriched matter back to the 
white dwarf via the accretion disk or as a descendant of a precursor binary that 
survived thermal timescale mass transfer. The scenario involving pollution of the 
secondary by past novae may be supported by the possible presence of a nova 
remnant-like nebula around EY Cyg. 

Subject Headings: Stars: cataclysmic variables, white dwarfs. Physical Processes: 
accretion, accretion disks 
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1. Introduction 

Dwarf novae are a subset of cataclysmic variables consisting of a white dwarf accreting 
matter from a Roche lobe-filling, low-mass secondary star. They undergo outbursts 
(brightness increases of 3-7 magnitudes) lasting for days to weeks separated by intervals of 
quiescence lasting weeks to years. The outbursts are thought to be caused by the release 
of gravitational potential energy when matter in the disk periodically accretes onto the 
white dwarf. The trigger for the accretion episode is widely held to be the disk instability 
mechanism (e.g., Cannizzo 1998 and references therein). During outburst, these systems 
are dominated by a luminous accretion disk, a high rate of accretion, and wind outflow. 
The spectra during quiescence are more difficult to interpret and generally exhibit a mix of 
absorption and emission features with a continuum energy distribution which may or may 
not be consistent with an accretion disk. 

In this paper, we utilize composite white dwarf photosphere plus accretion disk 
synthetic spectra to model the poorly understood far UV spectra of dwarf novae in 
quiescence. Our approach is described in Winter & Sion (2003) and was also applied to 
the HST spectra of RX And (Sion et al. 2001; Sepinsky et al. 2002). Our approach is 
most useful for dwarf novae in which we cannot easily identify the dominant light source 
during quiescence in the far UV. For the vast majority of dwarf novae, we do not know the 
accretion rate during quiescence, nor is the spectrum of the underlying white dwarf easily 
identifiable. Among the key questions we seek to answer are: What is the accretion rate 
during quiescence? Is this rate of accretion consistent with the disk instability model? How 
much flux does the white dwarf contribute to the FUV? How much flux does the accretion 
disk contribute to the FUV? How hot is the white dwarf? 

With these goals in mind, we expand our investigation of systems in quiescence, with 
the selection of the U Gem-type dwarf novae, EY Cygni. It exhibits a composition puzzle 
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shared by only a few other non-magnetic and magnetic CVs. It is well-known that the most 
prominent and common emission line in the far UV spectra of dwarf novae is almost always 
Civ (1548 A, 1550 A). EY Cyg appears to be an exception to this rule because it has very 
strong Nv emission and weak or absent CiV, as first reported by Winter & Sion (2001). 
Further evidence of this was seen in a recent study conducted by Gansicke et al. (2003), 
where examinations of HST/STIS spectroscopic snapshots of both BZ UMa and EY Cyg 
revealed both systems to have among the largest Nv/Civ ratios observed in CV systems. 
Could the depleted carbon be pointing to a different thermonuclear history for these dwarf 
novae, related either to donor star contamination by nova ejecta from the accreting white 
dwarf (see for example Sion et al. 1998), or nucleosynthetic processing in a donor secondary 
originally more massive than the accreting white dwarf that underwent a short phase of 
thermal timescale mass transfer (Schenker et al.2002)? Is there any property of these two 
systems which points to the reason for the peculiar Nv/CiV emission line flux ratio? 

The orbital period and recurrence time between outbursts of EY Cygni has recently 
been re-determined (Echeverria et al. 2002; Costero et al. 1998). The new orbital period is 
0.45932 days (compared with the previous value 0.2185 days reported by Sarna et al. 1995). 
The recurrence time of its outbursts, based upon long term AAVSO archival hght curve 
data, is 2000 days. Its light curve is highly variable and very unusual in that it appears to 
be shifted by 0.5 phase relative to usual dwarf nova light curves with similar parameters. 
All of this prompted Tovmassian et al. (2002) to obtain direct images of EY Cyg in Ha to 
search for a nebulosity or common envelope. Prom the radial velocity study (Echeverria et 
al. 2002) and the light curve estimate (Tovmassian et al. 2002, Costero et al. 1998, Smith 
et al. 1997), the inclination angle was found to be very low. They constrain i to be in the 
range of 13 to 17 degrees. 

An integral part of our synthetic spectral analysis is to derive a distance fom the model 
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fitting and compare it to estimated or derived distances to these systems as a consistency 
constraint on the goodness of fit. Distances can be derived from the absolute magnitude 
at outburst {My(^^ax)) versus orbital period relation of Warner (1995). This yields 505 pc 
for EY Cyg. Recent FGS parallaxes of dwarf novae by Johnson et al. (2003) have yielded 
a new relation between the My(^jnax) and orbital period, viz., My(^rnax) — 5.92 - 0.383 Porb 
(hrs). However, for EY Cyg, this relation yields 871 pc, which may be spurious due to EY 
Cyg's period being far above the upper boundary of the period gap. A reasonable upper 
hmit (700 pc) to the distance of EY Cyg is imposed by the fact that it is not embedded in 
the Cygnus Superbubble (see below) while a reasonable lower hmit (250-300 pc) is imposed 
from the spectral type of the secondary star. For EY Cyg, Tovmassian et al. (2002) find 
the distance to he between 400 and 500 parsecs. We adopt the mean of this range (450 pc) 
as the distance to EY Cyg. 

The observed characteristics of EY Cygni are given in Table 1 where, by column, 
we list: (1) orbital period in days, (2) recurrence time between outbursts in days (3) 
orbital inclination, (4) spectral type of the secondary, (5) mass of the primary in solar 
masses (6) mass of the secondary in solar masses, (7) apparent magnitude in outburst, and 
(8) apparent magnitude in quiescence. The references to the entries are listed below the table. 



2. FUSE observation of EY Cygni in Quiescence 

We obtained a FUSE spectrum of EY Cygni on 16 July, 2003 during its quiescence. 
The exposure time was 18,514 seconds through the LWRS aperture. The weakness of C III 
(1175) is consistent with the C depletion evident in the STIS spectrum (see below) and the 
lUE archival spectrum. 
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We also obtained an HST spectrum under the Cycle 11 snapshot program led by one 
of us (BTG). The details of the spectrum have already been published and can be found in 
Gansicke et al. (2003). 

3. Ground-based direct imaging 

The ground-based observations of EY Cyg were done at the 3.5 m WIYN^ telescope 
at Kitt Peak. The imaging camera with 2Kx2K CCD and a set of nebular filters was 
used. The observations in the filter and corresponding narrow band continuum were 
obtained on September 12, 2001. The other set of observations, in filters Sll and Olll, 
were obtained one year later on September 3, 2002 on the same telescope with the same 
instrument. Images in each filter were obtained in multiple short exposures in order to 
avoid saturation by bright stars. This method of combining multiple images also helped to 
get rid of cosmic rays in the best manner. The seeing was around 0.6 arcsec on average 
during both sets of observations. The seeing in continuum images was slightly worse than in 
Ha. In order to subtract continuum images from the emission line, we had to do a Gaussian 
filtering to degrade the Hq images to the level of continuum (about 0.1 arcsec). For the 
S II images we used the scaled background image of the previous year, since no continuum 
image was available in 2002. The table of optical observations is presented in Table 2. The 
data reduction of mosaic images, including bias processing, flatficlding, combining multiple 
exposures, and continuum subtraction, was done using the MIDAS astronomical package. 

The surroundings of EY Cyg show complex Ha clouds common to Cygnus. The images 
of EY Cyg in the Ha filter are presented in Fig 1. Fig la presents a larger area around the 

^The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indi- 
ana University, Yale University, and the National Optical Astronomy Observatory. 



Table 1. EY Cygni Parameters 



Porb tree i Sp. type(M2) Mi M2 Vmax Vmi, 



0.45932 2000 13-17 K7V 1.26 0.59 11.4 15.5 



Table 2. Log of the Optical Observations of EY Cyg 



UT Date 


Band 


Filter # 


^exp SGC. 


Number of exp. 


Sept 12 2001 


Ha 


W015 


300 


5 


Sept 12 2001 


Cont. 


KP1494 


300 


2 


Sept 03 2002 


Sll 


W017 


600 


3 


Sept 03 2002 


Oiii 


KP1467 


600 


2 
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object without continuum substraction. This panel also shows the immediate surroundings 
of EY Cyg after continuum substraction. The star is situated at the rim of a very circular 
ring of Hq! emission. We checked the possibility of EY Cyg having a proper motion that 
could move it to the center of that ring in a reasonable time stretch, but the position of the 
star does not change within the limits of the measurements in the past 50 years. It is very 
unlikely then that this ring is associated with our object. The other suspicious structure 
detected is in the form of a faint arc in the immediate surroundings of EY Cyg, shown in 
Fig lb. It is not homogeneous. The size is about 25 arcsec. The lower limit of the distance 
to EY Cyg from the secondary spectral class is at least 250 pc. The upper hmit is probably 
700 pc, otherwise the object would be embedded in the background Cygnus superbubble 
and could hardly sustain its X-ray luminosity (Bochkarev & Sitnik 1985). Therefore, a 
reasonable distance estimate to EY Cyg would be around 400-500 pc; this is also supported 
by the spectral class of the secondary. In this case, the arc would be comparable in size to 
the known nova shells of DQ Her and/or GK Per. The morphology of this formation is a bit 
different from larger scale strips crossing the region. Moreover, EY Cyg is situated in the 
center of the arc, suggesting that it may be associated with the object. The images in Sll 
show similar but much fainter structure around the star. We were not able to detect any 
emission in O ill from EY Cyg or from the larger area. The arc in S ll is only marginally 
seen, but its form coincides with that of the arc in Ha, which confirms its reality and 
suggests a possible association with EY Cyg. However, we cannot as yet definitively prove 
the association of the arc with the star. 

4. Synthetic Spectral Fitting 

Model spectra with solar abundances were created for high gravity stellar atmospheres 
using TLUSTY (Hubeny 1988) and SYNSPEC (Hubeny & Lanz 1995). We adopted model 
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accretion disks from the optically thick disk model grid of Wade & Hubeny (1998). Using 
lUEFIT, a minimization routine, both values and a scale factor S were computed 
for each model. The scale factor, normalized to a kiloparsec, can be related to the white 
dwarf radius through: Fx^obs) — ^T^iR^ / <P')Hx{rnodei), where d is the distance to the source, 
H is the Eddington (surface) flux and S — At:{R^ /cP). For each system, the best- fitting 
accretion disk model was combined with the best-fitting photosphere model to determine 
the contribution of both the accretion disk and the white dwarf. For this combined white 
dwarf plus disk model fitting, we took the best-fitting white dwarf model and combined 
it with the best-fitting disk model, but, in effect, modulated the accretion rate of the 
best-fitting disk model by multiplying the disk fluxes by a small numerical factor in the 
range 0.1 to 10.0. In this relatively narrow range, the scahng of the disk fluxes is roughly 
linear so the scaling approximates varying the accretion rate by a small increment and 
examining how it affects the value. We then found the best-fitting composite white 
dwarf plus disk model based upon the minimum value achieved and consistency of the 
scale factor-derived distance with the adopted distance for EY Cyg. 

The overlap in wavelength between the FUSE and STIS spectrum is 1150A to 1182A. 
The flux levels between the two spectra are mismatched by a factor of 1.85 with the FUSE 
spectrum having a higher flux level than the STIS spectrum. In order to get the best match 
of the spectra in the overlap region, we raised the entire STIS spectrum by a factor of 1.85, 
but ehminated the overlap region in the STIS spectrum and instead used the FUSE overlap 
fluxes in that range since they are more reliable there. While we cannot rule out that the 
mismatch is due to a STIS calibration problem, we believe it is more likely that EY Cyg 
was in a lower accretion state when the STIS spectrum was taken. Unfortunately, we have 
no optical observations which document brightness differences between the times of the 
STIS and FUSE spectra. If we had divided the FUSE fluxes by 1.85 to match the lower 
STIS fluxes, then a lower temperature from the overall fit would have resulted. Hence, the 
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temperature we have derived from raising all of the STIS fluxes to match the FUSE fluxes 
in the overlap range can be regarded as an upper limit. 

For EY Cygni, we constructed a grid of white dwarf models with fixed surface gravity 
log ^ = 9, Te///1000 = 18, 19,20, ....28, Si and C abundance 0.1, 0.2, 0.5, 1.0, 2.0 5.0 and 
Vrot = 100, 200, 300, 400, 600, 800. For the all of the accretion disk fitting, we adopted 

= I.2M0 and i = 18 degrees in order to be consistent with the published values and 
with the white dwarf gravity in the photosphere fits. For the accretion belt fitting, taking 
sin 18 ~ 0.31, Vheitsini = 2000 km/s. This high velocity region is rotating near the Keplerian 
velocity, much faster than the white dwarf rotation velocity of 400 km/s. However, the high 
velocity region covers a relatively smaller fraction of the white dwarf surface area and may, 
in reality, lie above the surface (e.g. a hot rotating ring) rather than being right at the 
photosphere. For the belt, we took log g = 6 and solar abundances for all elements. We 
tried a range of accretion belt temperatures T^eit/lOOO = 30, 31, 32, ....50 in steps of 1. 

The results of our fitting of a single temperature white dwarf to EY Cyg are somwhat 
encouraging in that the best-fitting model (i.e., least — 1-233) yielded a lower distance 
(301 pc) than our adopted distance of 450 pc but not unreasonably lower. We find that 
the if the FUSE +ST1S spectrum is due to the white dwarf alone, then Teff = 24, OOOK, 
log g — 9.0, V sin i — 600 km/s with Si abundance = 0.1 x solar and C — 0.2 x solar. Even 
though the fitting yielded a rotational velocity, this value must be viewed with caution due 
to the lack of well-defined line profiles in the spectra. Since the FUSE and STIS spectra 
were mismatched by a factor of 1.85 in the overlap region, the temperature we have derived 
should be considered an upper limit (see above). This best-fitting model is displayed in 
figure 2. Table 3 gives the fitting parameters, where, by column, we list (1) is the model 
type, (2) is the effective temperature, (3) the stellar surface gravity, (4) the rotational 
velocity in km/s, (5) silicon abundance, Si/H, relative to solar, (6) carbon abundance, C/H, 



relative to solar (7) the reduced value, (8) the scale factor, (9) the distance in parsecs. 
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Results 
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WD 


24,000K 


9.0 
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3.9 X 10-^ 
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On the other hand, wc tried accretion disk model fits alone. The best fitting disk 
model with M^j^ — 1.2Mq and i = 18 degrees was for an accretion a rate of lO'^^Mg/yr 
had a = 2.148, almost indistinguishable from the the white dwarf fit. However, the scale 
factor of the disk fit yielded a distance = 1.8 kpc! 

The disk model plus white dwarf giving the minimum x^ value was exactly the 
accretion rate derived from the lone accretion disk fit. In the combined fit, the accretion 
disk accounted for 96% of the FUV flux while the white dwarf accounted for 8% of the FUV 
flux. However, the distance derived from the scale factor of the combined flt was 1,060 pc, 
substantially larger than than our adopted most probable distance of 450 pc. 

Finally, we attempted flts involving two-temperature white dwarfs, in which a hotter, 
rapidly spinning region of small surface area (e.g., an accretion belt) is combined with a 
white dwarf rotating more slowly with a lower temperature. We also constructed accretion 
belt models with log c/ = 6, = 3400 km/s, T^eit/lOOO = 30,31, ...50. The best-fltting 
white dwarf plus belt model yielded the lowest x^ value (2.03) of all of the above single and 
two-component fltting experiments and a distance = 351 pc, within a 100 pc of the adopted 
distance of 450 pc. The white dwarf contribution to FUV flux is 37% while the accretion 
belt contributes 63% of the FUV flux. The accretion belt covers 27% of the white dwarf 
surface area, has a belt temperature of 36,000K, log g — & and Ybeit sin i — 2000 km/s. The 
best-fitting accretion belt plus white dwarf model is displayed in figure 3 and the fitting 
parameters are given in Table 4 where we list by column: (1) the model type (2) the white 
dwarf mass (3) the WD Tg// (4) the accretion belt temperature (5) accretion belt velocity 
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(6) percentage flux contribution of the WD (7) percentage flux contribution of the belt (8) 
the reduced x value (9) scale factor and (10) the distance in parsecs. 



Table 4: EY Cygni White Dwarf plus Accretion Belt Fitting Results 



Model 


Mtud Tg// 




^belt 


Flux (WD) 


Flux (Belt) S d 


WD+Belt 


1.2 22,000K 


36,000K 


2000 


37% 


63% 2.033 2.88 x 10"^ 351 



5. Conclusions 

Based upon the fltting experiments described in the preceding section, it seems clear 
that the FUSE plus STIS spectrum of EY Cygni in quiescence is best-represented by 
either a single temperature white dwarf with Te// = 24, OOOK or by a white dwarf with 
Te// = 22, OOOK and an accretion belt having Tij^u — 36, OOOK and Yheit sin i — 2000 km/s. 
Since the latter flt has a lower and yields a distance closer to 450pc, it is favored over 
the single temperature white dwarf fit. The accretion belt is the dominant source of FUV 
flux from the FUSE range through the STIS range with the cooler white dwarf latitudes 
contributing less than 40% of the FUV flux. 

In EY Cygni, which lies above the period gap, the T^s value appears to be lower than 
the typical 29,000 to 35, OOOK of other U Gem and Z Cam systems above the gap. Could 
there be any connection between the high mass white dwarf, lower Te// of EY Cygni and 
the abundance anomaly seen in the far-UV spectra? With a higher white dwarf mass and 
lower Teflf value, it is possible that the abundance anomaly (excess N, depleted C) may be 
associated with a greater system age (longer cooling time of the white dwarf) than other 
dwarf novae. If indeed the white dwarf core has cooled for a longer timescale, then the 
equilibrium Teff value in response to accretion is expected to be lower. This might suggest 
that an older system age and massive white dwarf allows for more classical novae episodes. 
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hence more CNO nucleosynthetic processing, consequently more contamination of the 
secondary's atmosphere during the brief common envelope stage of each nova. Such systems 
may be more likely to exhibit the apparent excess N and depleted C associated with CNO 
processing. Our direct images with a narrow Ha filter of EY Cyg and its surroundings show 
the nova- like appearing envelope around it, but were inconclusive as to whether the nebular 
emission around the object is associated with it. 

An alternative, more promising picture has emerged however. Recent work by Gansicke 
et al. (2003) suggest that the N/C anomaly seen in the dwarf novae BZ UMa, EY Cyg, 
IRXS J232953.9+062814, and now CH UMa (see Dulude & Sion 2002), may have its origin 
in a CV with an originally more massive donor star (M2 > l.SM©) which survived thermal 
time scale mass transfer (Schenker et al. 2002 and references therein). In such a system, 
the white dwarf would be accreting from the peeled away CNO-processed core stripped of 
its outer layers during the thermal timescale mass transfer. 
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Fig.la- A larger area than on the right is presented without continuum subtraction, 
showing the surroundings of EY Cyg. CV position represented by X. 

Fig. lb- A blow-up of the continuum subtracted image of EY Cyg. A faint arc crossing 
the larger scale strips of nebular emission can be seen to the north-east of the star. 

Fig. 2- The best-fit single temperature white dwarf model to the combined FUSE 
-I- HST/STIS spectrum of EY Cyg during quiescence. The white dwarf model has 

= 24, OOOK, log g = 8.96, rotational velocity Vsini = 600 km/s, with Si abundance = 
0.1 X solar and C = 0.2 x solar. 

Fig. 3- The best-fit combination of white dwarf plus accretion belt synthetic fluxes to 
the spectrum of EY Cyg during quiescence. The white dwarf model has Tgs — 22, OOOK, log 
g — 8.9 and the accretion belt has log g — 6, Vsini = 2000 km/s and covers 27% of the WD 
surface area. The top solid curve is the best-fitting combination, the dotted curve is the 
white dwarf spectrum alone and the dashed curve is the accretion belt synthetic spectrum 
alone. In this fit, the accretion belt contributes 63% of the far UV flux, and the white dwarf 
contributes 37% of the flux. 



